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In the vicinity of the QCD phase transition, critical fluctuations have been predicted to 
lead to non-statistical fluctuations of particle ratios, depending on the nature of the phase 
transition. Recent results of the NA49 energy scan program show a sharp maximum of the 
ratio of K + to 7r + yields in central Pb+Pb collisions at beam energies of 20-30 AGeV . This 
observation has been interpreted as an indication of a phase transition at low SPS energies. 
We present first results on event-by-event fluctuations of the kaon to pion and proton to 
pion ratios at beam energies close to this maximum. 

Recently the energy scan program at the CERN SPS providing Pb+Pb collisions at beam energies of 20, 
30, 40, 80 and 158 AGeV has been completed. The NA49 collaboration has measured identified charged 
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particle spectra at all provided beam energies. One of the most striking observations in this data set is 
a sharp maximum of the ratio of average K + to tt + yields in the range of 20-30 AGeV incident beam 
energy [1]. Within a statistical model this maximum can be interpreted as an indication of the onset of a 
deconfinemcnt phase transition. A potential test for the hypothesis of observing direct effects of a phase 
transition is the analysis of event-by-event fluctuations in the hadro-chemical composition of the particle 
source. It has been suggested that, depending on the nature and order of the phase transition, anomalies in 
the energy dependence of cvent-by-event fluctuations might be observed and they may allow to characterize 
the transition process [2-5] . 
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FIG. 1. Phase-space distribution of accepted charged particles at 20 AGeV (upper panel) and 158 AGeV (lower 
panel) . From left to right pions, kaons and protons are shown. No acceptance correction was applied to the particle 
ratios reported here. 



A previous analysis of the event wise ratio of charged kaons to charged pions ( [K + + K~]/ [tt + + ir~ ] ) was 
performed in Pb+Pb collisions at 158 AGeV. The non-statistical fluctuations of 2.8% observed in the data 
were found to be significantly smaller than those expected for an independent superposition of nucleon- 
nuclcon collisions [6]. This supports the interpretation that at the top SPS energy each collision samples the 
same flavor ratios as described in a grand canonical ensemble, combined with a smooth transition from a 
possible partonic state to the final state hadronic particle composition. The minimal fluctuations expected 
due to production of the final state hadrons via resonances completely exhaust the observed fluctuation 
signal [7]. 

In this article the extension of this analysis to lower incident energies is presented. In addition to the 
estimate of the [K + + K~]/[ir + + ir~] ratio also the ratio of [p + p]/[7r + + 7r~] is determined on an event-by- 
event basis. The analysis method is described in [6,8,9]. The extension of this method to fitting multiple 
particle ratios for each event and the robustness of the event-by-event ratio estimators against effects of 
very low particle multiplicities have been validated by extensive systematic checks and Monte Carlo studies. 

A detailed description of the NA49 experiment can be found in [10]. At the five available beam energies 
the 3.5% most central Pb+Pb collisions were selected based on projectile spectator energy. Figure 1 shows 
the phase space distribution of accepted pions, kaons and protons at 158 AGeV compared to the 20 AGeV 
data set. The shifting cut off for pions is due to the requirement of a minimum total momentum of 3 GeV/c 
in order to facilitate the particle identification. The relative width a, defined as a = RMS /Mean* 100 [%], 
of the measured event-by-event particle ratio distributions can be decomposed into three contributions: 

1. Due to the finite number of particles produced and observed per event, the ratio of particle mul- 
tiplicities measured event-by-event will exhibit statistical fluctuations with a width dictated by the 
individual particle multiplicities. 

2. Due to non- ideal particle identification these statistical fluctuations will be smeared by the experi- 
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mental dE/dx resolution and the event-by-event fitting procedure. 

3. Superimposed on the background of statistical and experimental fluctuations we may observe genuine 
non-statistical fluctuations. 

An accurate estimate of the contribution due to finite number fluctuations in the particle multiplicities and 
effects of detector resolution is obtained using a mixed event technique. Mixed events are constructed by 
combining particles randomly selected from different events, reproducing the multiplicity distribution of 
the real events. 
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FIG. 2. Distributions of the event-by-event [K + + K']/[tt + + tt 
(histogram). 
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FIG. 3. Distributions of the event-by-event [p+p]/[-K + ratio for data (points) and mixed events (histogram). 

These mixed events are then subjected to the same fit procedure as the real events. By construction, the 
mixed events on average have the same particle ratios as the real events, but no internal correlations. To 
quantify the relative contribution of dynamical fluctuations adyn, the r.m.s width a m i X of the mixed event 
distribution is subtracted from the r.m.s width data of the real event particle ratio distribution: 



a dyn = sign(crj ata - <J 2 mixed ) y/ >La ~ a Lxed\ 0) 

The systematic error of the fluctuation measurement is estimated by comparing results when using two 
different track quality cuts. The values of the fluctuation signal presented here are calculated as the 
arithmetic mean of both results. For illustration the distributions of the two studied particle ratios at 20, 
40 and 158 ^4GeV are shown in Figs. 2 and 3. The contribution from dynamical fluctuations obtained from 
Eq. 1 is plotted in the Fig. 4. Fluctuations of the K/tt ratio are positive and decrease with beam energy. 
In case of the [p + p]/ir ratio the width of the data distribution is smaller than the width of the distribution 
of mixed events. The dynamical fluctuations are negative. A negative fluctuation signal of the event wise 
[p+p]/n ratio can be understood, if resonance decays into pions and protons are considered. The magnitude 
of the negative fluctuation signal in the [p + p\/i: channel may be related to the relative contribution of 
resonance decay products in the final state of the collision. In order to estimate the significance of the 
observed fluctuation signals of the two ratios considered, we compare the data to a string-hadronic cascade 
model UrQMD [11]. In this model, by construction, no fluctuations due to a potential phase transition 
are present, while resonance decays are included as well as effects of correlated particle production due 
to quantum number and energy/momentum conservation laws. For this study, large samples of UrQMD 
events were generated at all five beam energies and then subjected to an acceptance filter modeling the 
NA49 detector system. The accepted final state particles were counted and the corresponding ratios were 
formed. 
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FIG. 4. Energy dependence of the event-by-event fluctuation signal of the [K + + K ]/[7r + + 7r ] ratio (left panel) 
and the [p + p]/[7r + + tv~] ratio (right panel). The systematic errors of the measurements are shown as gray bands. 



The energy dependence of the event-by-event [p + p]/ir ratio in UrQMD closely matches the energy de- 
pendence observed in the data, as shown in Fig. 4. This lends further support to interpreting the negative 
fluctuation signal as resulting from resonance decays. In case of the fluctuations of the event wise K/ir 
ratio, the energy dependence of the signal cannot be reproduced in the cascade model. UrQMD gives an 
energy independent fluctuation signal. Since the relative contribution of resonances changes dramatically 
with incident beam energy, we conclude that in the K/tt ratio resonances do not give a significant con- 
tribution to the fluctuation signal. The finite fluctuation signal in the UrQMD model can be attributed 
to correlated particle production due to conservation laws. In the data we observe a significantly smaller 
fluctuation signal at highest beam energies than in a cascade model. At 158 AGeV the fluctuation signal 
is consistent with calculations performed assuming a grand canonical ensemble without enforcing local 
conservation laws. Towards lower beam energies we see a steep increase of the fluctuation signal as shown 
in the left panel of Fig. 4. The increase of the signal goes significantly beyond the value seen in a hadronic 
cascade model, indicating the onset of a new source of fluctuations. Further theoretical calculations will 
be needed to evaluate the relevance of this observation for a possible interpretation in the context of a 
deconfmcmcnt phase transition. 
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